Big Bang Nucleosynthesis with an Evolving Radion in the Brane World Scenario 
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We consider the big bang nucleosynthesis (BBN) in the Brane world scenario, where all matter 
fields are confined on our 3-brane and the radion of the Brane evolves cosmologically. In the Einstein 
frame fundamental fermion masses vary and the results of standard BBN (SBBN) are modified. We 
can thus use the observational primordial element abundances to impose constraints on the possible 
variations of the radion. The possibility of using the evolving radion to resolve the discrepancies 
between the Wilkinson Microwave Anisotropy Probe (WMAP) and SBBN values of the baryon-to- 
photon ratio (r;) is also discussed. The results and constraints presented here are applicable to other 
models in which fundamental fermion masses vary. 
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I. INTRODUCTION 

The Brane world scenarios have been attracting much 
attention during recent years, being motivated by the de- 
velopments in superstring theory and acting as alterna- 
tives to supersymmetric or technicolor models in address- 
ing the gauge hierarchy problem. In such scenarios the 
matter fields are confined in our world, which is believed 
to be a 3-dimensional brane embedded in a 1 -I- 3 -I- n- 
dimensional bulk spacetime, while gravitation fills all the 
dimensions. In Ref. T| Arkani-Hamed, Dimopoulos and 
Dvali (ADD) proposed a model in which the fundamen- 
tal energy scale can be as low as TeV while the weakness 
of gravity is due to the large compactification radius of 
extra dimensions (millimeter in case of two extra dimen- 
sions, for example). A different setup was introduced by 
Randall and Sundrum (RS) later, in which our brane 
is one (the negative tension one) of the two boundaries of 
a 5-dimensional anti-de Sitter (AdS5) spacetime and the 
gravitation is localized around the positive tension brane 
while observers on our brane only measure a weak tail of 
it. In some works, e.g., H,^0, the various matter fields 
are also allowed to propagate along the extra dimensions 
and interesting consequences are discovered such as a dif- 
ferent unification scale and the localization of zero-mode 
fermions onto our brane. 

Brane model cosmology has also been investigated ex- 
tensively (see, ej., 0,01 for an introduction) . In particu- 
lar, it is found [g that the gravitational equations on our 
3-brane (in a 5-dimcnsional bulk, for example) appear in 
the form 



(4) 



TT^jy - E^vi (1) 



where t^^, is the energy momentum tensor confined on the 
brane and Mpj the fundamental (5-dimensional) Planck 
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mass; A4 is the effective 4-dimensional cosmological con- 
stant which is related to the bulk cosmological constant 
A and the brane tension A by 
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The 4-dimensional Newtonian constant Gn is defined by 



Gn 



47rA 



(3) 
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The TT^^ term in Eq. (I) is a quadratic function of t^^ 
and negligible when the energy density on the brane is 
much smaller than the brane tension, i.e., t^^/\ <g; 1, 
and the last term is a part of the high dimensional Weyl 
tensor carrying information of the gravitation outside the 
brane, which acts in 4-dimensional theory as a relativis- 
tic energy component — the so-called dark radiation. If 
the non-conventional field equations hold down to low en- 
ergies and the dark radiation term is nonzero, then the 
predictions of the standard cosmological model, such as 
BBN and CMB (Cosmic Microwave Background), might 
be modified dramatically. Consequently the observations 
of primordial element abundances and CMB power spec- 
trum would place stringent constraints on the parameters 
in Eq. (1) 

In Ref. the authors considered the effects of a non- 
standard cosmic expansion rate on the outputs of BBN; 
yet there is another source for modifying the SBBN yields 
in the Brane World scenario - namely the variation(s) of 
the fundamental constant(s). This is because, although 
in most extra dimensional models the moduli fields (for 
instance the radion in the brane model) are believed to 
have been stabilized before the commencement of the pri- 
mordial nucleosynthesis, it is not necessary to dismiss the 
possibility of cosmological evolution of the moduli fields 
in order to avoid the long range forces and deviations 
from general relativity (GR) |ll|. It is conceivable that 
there is a cosmological attractor mechanism as in some 
scalar-tensor theories, driving the model towards the con- 
ventional GR [l^]; it is also probable that the radion it- 
self is a chameleon field acquiring a large effective 
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mass via its self-interaction together with the interac- 
tion with matter and thus evading the constraints from 
local gravitational measurements while evolving cosmo- 
logically 0, Such a time evolution of the radion 

field will lead to variations of fundamental constant (s). 
BBN turns out to impose the most stringent constraints 
on the possible time evolution of the extra dimensions be- 
cause it is more sensitive than, say, CMB to the changes 
of the fundamental constants. Therefore, it is our in- 
terest in the present work to find out how the outputs 
of BBN are infiuenced by an evolving radion field (See 
[TgL \vA ITsI for the BBN constraints on other specific 
models and 0, |23, 0, 113, 113 etc. for constraints on 
the variations of individual fundamental constants such 
as the fine structure constant and the strong coupling). 

This article is arranged as following: in Sec. II we dis- 
cuss how an evolving radion field changes the fundamen- 
tal constant (s) and thus leads to modified predictions of 
BBN. We will work in the Einstein frame where Gn stays 
unaltered and show that the only one varying fundamen- 
tal constant in the model is the Higgs vacuum expecta- 
tion value (VEV) v = {H). The BBN constraints on a 
changing Higgs VEV have been investigated previously 
in |24l l25j | and more recently in , but ours differs from 
the previous works in several ways: firstly, we consider 
some effects ignored in previous works, such as the ra- 
diative and Coulomb corrections to the neutron lifetime 
(and weak interaction rates), their radion dependence 
and the radion dependence of the p(n, D)^ cross section; 
secondly, we have done a complete likelihood analysis us- 
ing the recent compilation on the nuclear rates and un- 
certainties ll^, 113 and measurements of the primordial 
abundances of D, ''He and ^Li; thirdly, we show that the 
^Li yields could be changed significantly with a small evo- 
lution of the radion, tending to reduce the inconsistency 
between the SBBN and WMAP-implied values of 77, the 
baryon-to-photon ratio 29] . Because of the accuracy of 
the WMAP result, we also use it to constrain the varia- 
tion of the radion. Our numerical results obtained from 
a modified BBN code is presented in Sec. HI and then 
Sec IV is devoted to a discussion. Throughout this work 
we will assume three species of massless neutrinos (or 
negligible neutrino masses) and zero chemical potentials 



for neutrinos as in SBBN. 



II. THE INFLUENCES OF AN EVOLVING 
RADION ON BBN 

In this section we discuss the influences of an evolving 
radion on the primordial nucleosynthesis. We first write 
down the low energy action in the 4-dimensional effec- 
tive theory and show that only the Higgs VEV (and thus 
fundamental fermion masses) changes. For simplicity we 
will concentrate on large flat extra dimensions; the more 
realistic models with warped extra dimensions are dis- 
cussed in the literature, e.o., [13. (In J 2] the authors 
considered a general class of two brane models, includ- 
ing the RS model as a special case; the time dependence 
of fermion masses in the Einstein frame and time inde- 
pendence of gauge couplings can also be found there.) 
Then we briefly discuss the implications on BBN. The 
derivations in this section closely follow our recent work 



A. The Low Energy Effective Action 

Let us start from a general 4 + n-dimensional model, 
with n being the number of the large extra dimensions. 
The full line element is given as: 

ds^ = GABdX^dX^ = g^^dx^'dx" + habdy^dy^ (4) 

where fj,, v — 0, 1, 2, 3 label the four ordinary dimen- 
sions, a, & = 4, • ■ •, 3 -|- n denote extra dimensions and 
A,B = 0, 1, 2, 3+n describe the whole spacetime. (We 
shall not consider cross terms such as Gafj. in Eq. (4)). 
The extra dimensions are assumed to compactify on an 
orbifold, and their coordinates ya take values in the range 
[0, 1]. The quantities hab have dimensions of [Length]^ 
since ya are dimensionless in our choice. 

Then the effective 4-dimensional action in the gravita- 
tional sector can be obtained by dimensionally reducing 
Eq. (4) as: 



ravity 



d^xd"y^\g\ 



|G|i?4+,i [G 

W\ 



Vn 



-h^-'d^Kb ■ h'^d^^Kd 



(5) 



in which \g\, \h\ and \G\ are respectively the determinants 
of the metrics of the ordinary dimensions, the extra di- 
mensions and the whole spacetime. Ri[g] and i?4_|_„[G] 
are the Ricci scalars of the ordinary 4 and the total A + n 



dimensional spacetimes. K4, K4+„ are related to the 4 and 
4 + n dimensional Planck masses through kI = 2M^i4 

and K4^„ = 2A/pj^4, while they themselves are con- 
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nected by a volume suppression 



4+n 



V, V being and choose the field i? to satisfy 



a measure of the extra space volume whose present-day 
value is denoted by Vq in Eq. (5) (Note that because 
of the specified choice of Vq and because the higher di- 
mensional quantity K4+„ is treated as a constant, the K4 
above also takes its currently measured value and is a 
constant rather than a variable). 

The effective Ricci curvature term is not canonical in 
Eq. (5); to make it so, let us take the conformal trans- 
formation 



(6) 



Vo 



(7) 



Then we obtain the effective 4-dimensional gravitational 
action in the Einstein frame: 



Sg 
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(8) 



r 



We shall make a further assumption that the extra di- 
mension(s) are homogeneous and isotropic, i.e., the met- 
ric of the extra space takes the following form: 



Next we turn to the matter sector of the effective ac- 
tion, firstly for the scalar fields. The action of a brane 
scalar field ^ is given by: 



hab = diag(-6^, • • •, -6^), 
and then the action Eq. (8) could be rewritten as 

'S'Gravity = j d'^X^/\g 

by defining a new scalar field, the radion a: 



+ lg'"'d^ad,a 
K4 2 



1 n + 6" 

K4\ 11 Vo 



(9) 



(10) 



(11) 



(12) 



because the n large extra dimensions are inscient to the 
field and have been integrated out. Then the same con- 
formal transformation given in Eqs. (6) and (7) trans- 
forms Eq. (12) into the following form: 



-I 



d'^x^/\g\ { -exp 



n + 2 



9^''d^(t>d^^ - exp 
I 



-2k 



n + 2 



u{4>) . 



(13) 



Note that from now on we will use k instead of K4 for where the potential U {ip) is given by 
simplicity. The kinetic term of the scalar field in Eq. (13) 
could be made canonical by redefining a new scalar field 

as " vi-/ 2""*" + 



U iv) = ^exp 



(16) 



ip = exp 



"2''Vn + 2' 



provided that the original potential tJ{<j)) takes the fol- 
(14) lowing form: 



and the action becomes (up to the higher order derivative 
term): 



= / 



d x^/\g\ 



(15) 



(17) 



The same technique could be applied to gauge fields, 
whose action before conformal transformation is given as 

^Gauge = - / d^X^l^F^^'^P;^,, (18) 
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where g is the gauge coupHng constant (a tilde is used 
to distinguish it from g, the determinant of the metric) 
and F^j^ are corresponding gauge field strengths. With 
the conformal transformation Eq. (6) the above action is 
changed into 

1 



^Gaugc = -J d'^x^lgl—F^'^'^F;^, (19) 

where g is unchanged because F^i, — have zero con- 
I 



formal weights and are unaltered under conformal trans- 
formations. 

We can also obtain the effective action for the Dirac 
fermion field in a similar way. Starting from the brane 
action: 

S** = y" d^x^\ {i*7'^D^* - m^^) (20) 
and using the conformal transformation Eq. (6) we get 



51 \ exp 



3k 
T 



z^'7''D^^ - exp 



'2k 



(21) 



r 



To make the kinetic part of the fermion action Eq. (21) 
canonical, we rescale the field as 



i]} = exp 



3k 



(22) 



and then using the conformality of the coupling of mass- 
less Weyl fermions, we rewrite Eq. (21) as: 



in which 



\g\ {i4>rD,,^j - m{a)i>^] , (23) 



r7i(cr) = exp 



(24) 



Note that our results above are equal to those of |30j 
when there are no universal extra dimensions (i.e., n = 
in their model). 

It is apparent from Eq. (24) that if the radion a evolves, 
then the fermion masses also vary. As in the standard 
model, we assume that the fundamental fermion masses 
are generated by the Higgs mechanism and that the Higgs 
potential takes the form as Eq. (17). We find from 
Eq. (16) that the Higgs VEV (w), which is obtained by 
simply minimizing the Higgs potential, has the same ra- 
dion dependence as the fermion masses (Eq. (24)) so that 
the Yukawa couplings are independent of radion. Recall- 
ing also that the gauge couplings would not depend on 
radion (Eq. (19)), which is because of, as pointed above, 
the conformal invariance of the gauge kinetic term and 



J 



this is true in general scalar tensor gravity theories [Tg . 
So here we indeed encounter a model in which the only 
varying fundamental constant is v. On the other hand, if 
the gauge fields exist in the full dimensions as suggested 
by the universal extra dimensional scenario then the 
above conclusion does not hold anymore and t he g auge 
couplings will be dependent on the radion field [l8j. 
B. BBN with a Varying Higgs VEV 

The primordial nucleosynthesis with a varying Higgs 
VEV has been discussed previously in 0, |23, llg , and 
here we shall only briefly identify the places where effects 
of Higgs VEV enter. We also omit the introduction to the 
BBN theory and observation and refer these to existing 
literatures, e.g., [H HI [H H HI S El . 

A different v from its standard value will modify BBN 
mainly in two aspects, the weak interactions and the nu- 
clear reactions. Consider firstly the weak interactions 

Ue + n ^ p + , 
+ n ^ p + i7f,, 
n ^ p + + i7e, 

which interconvert the neutrons and protons. Roughly 
speaking, they are important for BBN because they de- 
termine the neutron density at the beginning of BBN, 
and, because nearly all the neutrons are incorporated 
into ''He at last, they are crucial for the final ^He output 
(for more detailed discussions see e.g., H^l). The rates 
of these n ^ p interactions could be well described as 



r(n ->p) 



A 



de 



A 



de 



e(e-(z)2(e2_ 1)1/2 



[1 + exp(-eze)] {1 + exp [(e - q)z^]} 
e(e + g)2(e2_ 1)1/2 



[1 exp(eze)] {1 + exp [-(e -I- q)z^]}' 



(25) 
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T{p n) 



A 



A 



de- 



e(e-g)2(e2_ 1)1/2 



[1 + exp(eze)] {1 + exp [{q - e)z^]} 
I [1 + exp(-eze)]{l+ exp [(e + 5)^1/]}' 



de- 



(26) 



r 



using the Born approximation, where we have defined 



dimensionless quantities q 



o/rrie, e — Ec/rrie and 



Zv = nie/Ti, and Zg = m^/T with me being the electron 
mass, m„p the neutron-proton mass difference and T^, (T) 
the temperature of the neutrinos (the electromagnetic 
plasma). Here A is a normalization factor determined 
by the requirement that at zero temperature T{n — *■ p + 
e~ +z7g) = with T„ the neutron lifetime, which means 
that: 



A 



where 



de eie~qne'-l) 



a/2 



(27) 



(28) 



/^From Eqs. (25)-(28) we can see that v determines the 



weak rates through Gp, 



and 



respectively. The 



dependences of Gp and nie on v can be simply parame- 
terized as 



Gp oc — , 
rrie oc V. 



(29) 
(30) 



The neutron-proton mass difference m„p is a sum of the 
electromagnetic contribution (~ —0.76 MeV) and the u— 
d quark mass difference (- 2.053 MeV) [13 , for which the 
former is unchanged in the present model because the 
(electromagnetic and strong) gauge couplings are radion 
independent while the latter is proportional to the Higgs 
VEV. So iTLnp could be described as 



in which 



m„p = 2.053/9 MeV - 0.76 MeV, 



wbbn 

WNOW 



(31) 



(32) 



is the ratio between the Higgs VEVs at the BBN era and 
at present. TO„p is important also because it determines 
the equilibrium neutron-to-proton ratio through |3l| 



rin 

— = exp 



T 



(33) 



Although the Born approximation in Eqs. (25) and 
(26) captures the essential features, various corrections 
are needed for more accurate estimations (see e.g., 
[3^ m, E3)- -f'oi' example, the accuracy of the the- 
oretical value of Tn from Born approximation is ~ 7% 



while including the zero-temperature and Coulomb cor- 
rections reduces it to be ^ 0.1% These two correc- 
tions are the most important ones and the radiative cor- 
rection itself depends on m„p (and thus on p); in our cal- 
culation we numerically integrate Eqs. (25) and (26) and 
include these corrections explicitly. We have neglected 
other corrections for simplicity and as the accuracy is 
adequate for our purpose here (see [s^, IH, for 
discussions on this issue). 

Next we turn to the nuclear reaction sector. The 
change of v will modify the pion mass m^r, which is re- 
lated to the light quark mass iriq [4ll |: 



oc rriq oc 



(34) 



According to recent investigations of the dependence of 
the nuclear potential on the pion mass the 
deuteron binding energy 8^ shows a strong decrease 
when m,r increases; although the relationships between 
Bd and as obtained in |42l l43j| have large uncertain- 
ties, they could be well approximated by 



Bd.o 



(r + 1) 



(35) 



within the small range of (or p) we are working with 
[2^0, in which B^ q and m7r,o denote their present-day 
values and the central value of r ranges from 6 42] to 10 
|43l |. In the following we shall work with r = 10 and 6 in 
parallel and show that the resulting qualitative features 
are the same. 

The deuteron binding energy plays an important 
role in BBN because it determines whether a signifi- 
cant amount of D could be produced (via the reaction 
p{n,D)j), which then leads to the synthesis of heavier 
elements. Because of the huge number of photons (recall 
that 1] ^ 10^^°!), the inverse process 7(D,n)p could be 
active down to low energies (~ 0.08 MeV) and sensitive 
to the variations in Bd |2^ H^. Besides, the forward 
cross section of p(n,D)7 will also depend on Bd 0|i to 
estimate this dependence, we adopt the cross section for- 
mula calculated using the effective field theory without 
pion . The authors of E3| gave an expression 
of the p{n, 0)7 cross section in terms of several param- 
eters in which the deuteron binding energy Bd has by 
far the largest leverage. In our modified BBN code we 
average the cross section given in [Tg] over the thermal 
distribution of the particles numerically and obtain the 
reaction rate as a function of the temperature |48j |: we 
have checked that the results obtained in this way are 
essentially identical to those obtained using the fitted 
p(n,D)7 rates of [Ulil. 
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FIG. 1: The primordial abundances of *He, D and ^Li as a 
function of the baryon-to-photon ratio rj, for various values of 
p Eis indicated by the legend in the lower panel and r — 10. 
Also shown are the observational la ranges of the abundances 
of *He from [s^ (grey region) a nd [sgll (dark grey region) , D 
from l60j(grey region), ^Li from '6l1 (grey region), as well as 
the baryon-to-photon ratio implied by WMAP [29||. ryvvMAP = 
(6.14 ±0.25) X 10-^° (the vertical lines). 



For the other reactions which are of importance to 
the D yields, such as D(D, nfRe, D{D,p)^R, D(p, -ffRe, 
D(^H,n)''He and D('^He,p)^He, their cross sections may 
also be dependent on niT^ or but we have no similar 
effective field theory calculations on them; so we ignore 
them in the present article. If we simply take these cross 
sections to be related to the size of the deuteron radius, 
i.e. a (X 1/Bd, as in Ref. ^2^, then in the interested range 
of T] and p (See FIG. 5 below) we find modifications to the 
final outputs of D, "^He and "^Li by <- 6%, <- 7 x 10""* 
and <~ 3% respectively, all well lying within the corre- 
sponding ItT observational uncertainties (~ 16% for D, 
- 0.84% for ^He and - 21% for ^Li, see Eqs. (36)-(38) 
below). In the present work we do not adopt the 1/Bd 
parametrization of these cross sections due to the lack of 
more explicit expressions; rather we treat the above esti- 
mations as possible errors introduced by neglecting their 



FIG. 2: Same as Fig. 1, but here r — 6 and the values of p 
are different. 



m^r or Bd dependences and emphasize that further re- 
lated theoretical calculations are needed to reduce these 
errors. 

Note that the variation of the Higgs VEV probably will 
also cause modifications in the binding energies of other 
nuclei and thus change the corresponding reverse reaction 
rates. However there are again no explicit calculations of 
these modifications. Fortunately these effects are small 
because the abundances of these heavier nuclei fall far 
below their nuclear statistical equilibrium (NSE) values 
long before the beginning of BBN and the reverse reac- 
tions are effectively switched off. For example, according 
to our constraints on p, i.e., p ~ 1.01 (see below), the 
deuteron binding energy Bd will be decreased by ~ 5%. 
Let us suppose that the other binding energies have vari- 
ations of similar magnitude, since different nuclei may 
appear in the left-hand and right-hand sides of a reaction 
{e.g., ^He and '^H in •^He(n,p)'^H) and we do not know 
whether the variations of their binding energies correlate 
or anti-correlate. We thus arrived at a conservative es- 
timation of ^ 10% for the Q-values of these reactions. 
If the Q-values of all the reactions relevant to BBN are 
changed by 10%, then the abundances of D, '*He and ^Li 
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FIG. 3: Upper panel: individual contours of D (dotted), *He 
(dashed) and 'Li (dot-dashed) at 68% confidence level. Lower 
panel: joint contours of D+'^He+'^Li at 68% (dark grey re- 
gion) and 95% (grey region) confidence levels; the best-fitting 
parameters (r?,p) ~ (6.28 x fO"^", 1.009) are denoted by the 
white circle and the case of SBBN by the horizontal dash- 
dotted line. In both panels the range of r;wMAP is represented 
by vertical solid lines. Here r is taken to be 10. 



would be varied by no more than ^ 0.1%, ~ 5 x 10~^ and 
^ 0.5% respectively, in the range of rj of interest |50| |. 

There is also a modification to the cosmic expansion 
rate in the present model. As the Newtonian constant 
Gn is unchanged in the Einstein frame, this modification 
originates from the variation of energy density compared 
with SBBN — due to the modified electron and positron 
masses as described in _26J . We include this effect in 
the calculation although it is small compared with others. 



III. NUMERICAL RESULTS 

We have incorporated all the effects discussed in the 
previous section into the standard BBN code by Kawano 
|51| and used it to obtain the numerical results presented 
in this section. The nuclear reaction rates are taken from 
the NACRE compilati on 1271 (for the rates not given in 
|27j we adopt those of |32j|) and their uncertainties from 
the work of Cyburt, Fields and Olive 28] . For the present 



FIG. 4: Same as Fig. 3, but for r = 6. The best fitting point 
here is (7?,p) ~ (6.67 x 10"!°, 1.020). 



neutron lifetime we use the value suggested by the Par- 
ticle Data Group HI, r^'^ = 885.7 ± 0.8 s. 

In Fig. 1 we plot the abundances of the light nuclei 
"^He, D and ''Li with a modified Higgs VEV (v), which is 
characterized by p, in the case of r = 10. It is apparent 
that if V is larger at the BBN era (p > 1), then the deu- 
terium output increases while the '*He and ^Li (for ^Li 
we only consider the larger-ry-case implied by the WMAP 
result f29|) yields decrease compared with SBBN. The 
behavior of the ^He output is a consequence of several 
effects: firstly, a p-value larger than 1 leads to a larger 



(from Eq. (31)) and a smaller neutron density at the 
beginning of the nucleosynthesis (from Eq. (33)), thus fi- 
nally to a smaller '^He output; secondly, the weak inter- 
action rates will be smaller than those in SBBN, leading 
to increased final ^He abundance. Thirdly, the deuteron 
binding energy becomes smaller (from Eq. (35)) than 
in SBBN so that the nucleosynthesis will commence at a 
later time and become less efficient, producing less ^He 
and ^Li while leaving more D unprocessed (there is an 
extra decrease in the forward rate of p{n, D)"f, whose in- 
fluence is small compared with the effect of commencing 
BBN later). The sum of these effects is a smaller ''He 
abundance (see Fig. 1 of |26j | for a more explicit compar- 
ison among these effects). For comparison we also plot 
in Fig. 2 the case of r = 6, where the effect of Bd be- 
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FIG. 5: The likelihood functions as functions of p alone with 
WMAP constraint on rj imposed, for r = 10 (the solid curve) 
and r — 6 (the dot-dashed curve). Both curves are normalized 
such that the areas under them are unity. 



comes weaker than r = 10 but the essential features are 
the same. 

The results in Figs. 1 and 2 suggest that the discrep- 
ancy between the SBBN theory and WMAP observations 
is reduced if p is greater than 1 [s^l • It is well known that 
adopting the WMAP-implied value of 77 "2^, the SBBN 
can reproduce the observed D abundance, but more ''He 
and '^Li than their observational abundances. One possi- 
bility is that this discrepancy is due to some systematic 
errors in the *He and ^Li observations. However, it may 
also originate from some new physics, such as a varying 
fine structure constant [U, an altered deuteron bind- 
ing energy more than three relativistic (effective) 
neutrino species [H^l, a lepton asymmetry [s^, or some 
combination of them [s^ . etc. Our results show that we 
might add another candidate, a varying Higgs VEV, to 
this list. To see this point more quantitatively, and to de- 
rive a constraint on the possible evolution of the radion, 
next we will give the likelihood analysis of this model. 
For this purpose, we choose to use the linear propaga- 
tion approach proposed by Fiorcntini et al. [s^l and then 
generalized by Cuoco et al. |58| to estimate the error 
matrix. The observational abundances and uncertainties 
of ^He, D and ^Li are taken from Luridiana et al. [s^ . 
Kirkman et al. |60| and Bonifacio et al. respectively 
as: 



Y^" = 0.2391 ±0.0020; 



(D/H)°'^'' 
(^Li/H)°'''^ 



2.78t°i^ X 10- 



2.19: 



-0.46 
-0.38 



X 10 



-10 



(36) 
(37) 
(38) 



(Note that the ''He abundance is quantified by its mass 
fraction.) These observational results are also shown in 
Figs. 1 and 2 as grey regions, and in Fig. 1 we also show 
the earlier wider range of ^He given by Olive, Steigman 



\robs 
^ P 



0.238 ± 0.005, 



(39) 



by the dark grey region as a comparison, even though we 
will not use it in our likelihood analysis. 

The results for r = 10 are plotted in the upper panel 
of Fig. 3, where we treat 77 and p as two free param- 
eters and give the 68% C.L. contours from the abun- 
dances of D, '^He and ^Li individually. It is apparent 
that with the presence of an evolving radion field (or 
equivalently a varying Higgs VEV so that p 7^ 1) the ob- 
servational abundances of these 3 elements can be com- 
patible at the 1(7 level and consistent with the WMAP 
result. This qualitative conclusion is then confirmed by 
the lower panel of Fig. 3, in which the joint constraints 
from D+^He+^Li at 68% and 95% C.L. are shown. Our 
best-fitting point in the parameter space (the white cir- 
cle) is at p ~ 1.009 and 77 ~ 6.28 x 10"^°. This cen- 
tral value of 77 is slightly larger than that of 77WMAP 
(wmap = 6.14^0 25 X 10^^") but falls into the la range 
of the latter. In contrast, the SBBN (horizontal line) 
is only allowed marginally at 95% C.L. and should be 
further excluded at the same level with 77WMAP taken 
into account. A similar situation is presented for r = 6 
(Fig. 4), while the details are a little different: here the 
best-fitting point - (6.67 x 10'^°, 1.02) fies slightly out- 
side the la range of ?7vvmap and SBBN is excluded at the 
95% C.L. using the observational primordial abundances 
only. In any case we see an improvement compared with 
SBBN if p is appropriately chosen (note also that, if the 
1/Bd scaling of the cross sections discussed in Sec. II B 
is a good approximation, then p > 1 will also increases 
the rates of these reactions by reducing Bd, making the 
destructions of D more efficient; this will lead to less D 
yields and shift the D and joint contours leftward). 

In some cases we may be interested in how p alone is 
allowed to change, and for this purpose it proves to be 
convenient to use the constraint on 77 from the WMAP 
measurement (the 77WMAP given above) by virtue of its 
high accuracy. We shall construct a likelihood function 
of p only as [i^ : 



L{p) oc 



POO 








/ exp 

J —00 


[ J 


exp 


2 



(40) 

in which rj^ and a^ are respectively the central value and 
la range for rj given by WMAP; xiViP) is a function of 
both rj and p calculated using the same method as above 
(for D±''He±^Li). The results are shown in Fig. 5, where 
we give the cases for both r = 10 and r — Q and we 
have found that the 95% C. L. ranges for p in these two 
cases are p = 1.0089 ± 0.0012 and p = 1.0182 ± 0.0024, 
corresponding respectively to changes in the extra space 
volume V as: 

1.52 X 10^2 < — < 1.99 X 10^2 (41) 
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and 

3.09 X 10^2 < ^ <4.00 X 10"^ (42) 

where AV = Vq — Vbbn is the difference between the ex- 
tra space volumes at present and at the BBN era. The 
allowed variation of p (or equally the Higgs VEV) we find 
here is comparable to the value quoted in but the 
full likelihood analysis using the ttwmap constraint and 
the new observational abundances of D, "^He and ^Li pre- 
sented here strongly disfavor p < I. Furthermore, as seen 
from Fig. 5, the ranges for p with r = 10 and r = 6 do 
not agree with each other at the 95% C. L., even though 
in both cases the qualitative features are the same. This 
simply reflects the demand of a more precise understand- 
ing of how the deuteron binding energy depends on the 
light quark masses. 

If no constraint on 77 is used, we obtain looser con- 
straints on p, namely p = 1.0091 ±0.0071 for r = 10 and 
p = 1.0198 ± 0.0140 for r = 6 (again at the 95% C.L.), 
which are comparable to the results above. In this case 
a variation of the Higgs VEV as large as ^ 3.4% com- 
pared with its present value is still allowed at the BBN 
era (redshift z - 10^ - 10^°). 

IV. DISCUSSIONS AND CONCLUSIONS 

In summary, we have considered in this article the 
possible implications of a cosmologically evolving Brane 
moduli field on the outputs of the primordial nucleosyn- 
thesis. We begin with the discussions on how the funda- 
mental constants are modified if the radion field of the 
Brane varies and how these modifications would affect the 
BBN yields. Some effects not considered previously are 
included in this procedure. Then we present a likelihood 
analysis using the recent compilation of the various nu- 
clear reaction rates and observational abundances of 
the hght nuclei D/He and "^Li ,59,M Ji]- The WMAP 
constraint on ij j29l | is also adopted in determining the 
constraint on the quantity p, which characterizes the vari- 
ation of the Higgs VEV (and thus the evolution of the 
radion field in the present scenario). We find that the 
BBN yields could be changed in such a way that the dis- 
crepancy between SBBN and ttwmap might be reduced, 
provided that the Higgs VEV was slightly larger at the 
BBN era. This conclusion is robust within the present 



theoretical uncertainty of r (6 10), and the errors in- 
troduced by neglecting the p-dependences of other nu- 
clear reactions are estimated to well fall within the Ict D, 
*He and ^Li observational uncertainties. However, fur- 
ther developments in the understanding of these topics 
might help constrain p more accurately. The constraints 
we obtain in this work are also applicable to other mod- 
els in which the Higgs VEV (and fundamental fermion 
masses) varies cosmologically, and could be used to con- 
strain parameters in such models. 

It would be also interesting to consider other implica- 
tions of such an evolving radion field. One example is its 
influences on the CMB; this was discussed in Refs. p^l63 | 
and the authors found that these influences were mainly 
through the variation of the electron mass mg (see 62] 
for more details). This might provide another constraint 
on the radion fleld, but it is much looser than the one we 
obtain and at a different cosmic era (z - 10^). Another 
potential implication of interest is the possible existence 
of di-proton or di-neutron. This is because the bindings 
of the di-nucleon systems are mainly contributed by the 
long range pion force which depends on pion mass ttItt. 
If m-rr was smaller at the time of BBN, then the bind- 
ings of di-proton/di- neutron would be enhanced, and if 
the di-proton becomes bound, it will open a rapid chan- 
nel for the hydrogen fusion '63] and thus be catastrophic 
to the stellar lifetimes. The conditions for the di-proton 
and di-neutron binding were studied in [6^ l65l | . In the 
present work, we see that p > 1 is preferred by our BBN 
analysis, which means that TOtt.bbn > 'tIjt^now and the 
di-nucleon systems would be kept unbound during the 
BBN era. As the last example of possible implications, 
let us consider the stability of '^He which, if becoming 
stable, would fill the mass-5 gap in the BBN nuclear re- 
action chain, drastically enhancing the production of ^Li. 
In 23| the authors investigated this question and found 
that 5 (mq / Aqcb)/ (n^q/^QCi)) ^ —0.1 was required to 
prevent '"^He from becoming stable, which is obviously 
satisfied by our constrained results. 
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